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Abstract 

In dislocation-free martensites the components of the elastic strain tensor are constrained by 
the Saint- Venant compatibility condition which guarantees continuity of the body during external 
loading. However, in dislocated materials the plastic part of the distortion tensor introduces a dis- 
placement mismatch that is removed by elastic relaxation. The elastic strains are then no longer 
compatible in the sense of the Saint- Venant law and the ensuing incompatibility tensor is shown to 
be proportional to the gradients of the Nye dislocation density tensor. We demonstrate that the 
presence of this incompatibility gives rise to an additional long-range contribution in the inhomo- 
geneous part of the Landau energy functional and to the corresponding stress fields. Competition 
amongst the local and long-range interactions results in frustration in the evolving order parameter 
(elastic) texture. We show how the Peach-Koehler forces and stress fields for any distribution of 
dislocations in arbitrarily anisotropic media can be calculated and employed in a Fokker-Planck 
dynamics for the dislocation density. This approach represents a self-consistent scheme that yields 
the evolutions of both the order parameter field and the continuous dislocation density. We illus- 
trate our method by studying the effects of dislocations on microstructure, particularly twinned 
domain walls, in an Fe-Pd alloy undergoing a martensitic transformation. 

PACS numbers: 81.30.Kf, 63.70.+h, 61.72.Lk, 05.70.Fh, 05.10.Gg 
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I. INTRODUCTION 



A mesoscopic description (nano to micrometer) of physical processes in solids, where 
atomic length scales merge with those of the continuum, represents a crucial and perhaps 
most challenging aspect of understanding material behavior. This arises, for example, during 
displacive (martensitic) phase transformations where the distortions associated with the 
strains in unit cells and intra-unit cell displacements (or shuffles) propagate over larger 
distances so that competing long-range effects lead to the formation of inhomogeneities 
such as interfaces, spatially correlated domains and complex microstructure. It is the least 
understood regime compared to the atomic and continuum scales because the simplifications 
and advantages of theory in handling small/large length scales and fast/slow time scales no 
longer apply The predictions of atomistic models become typically invalid at length scales 
larger than about a nanometer, whereas the lack of the detailed description in continuum 
theories makes them inappropriate for studies of physical processes occurring at length scales 
below a millimeter. 

Our focus will be on materials undergoing structural phase transformations that represent 
an important and broad class of advanced materials frequently utilized in state-of-the-art ap- 
plications such as surgical tools, artificial muscle fibers, aerospace and robotic applications, 
and novel microelectronic devices. Examples include shape memory alloys (Ni-Ti, Fe-Pd, 
Au-Cd, Cu-Al-Ni), nuclear materials (Pu, U-Nb), ferroelectrics with spontaneous polariza- 
tion (BiFeOa, LuFe204, LiC^C^), strain-induced ferroelectric perovskites and high-k insula- 
tors (BaTi0 3 , SrTi0 3 , PbTi0 3 , LaA10 3 ), magnetic shape memory alloys (Fe-Pd, Ni 2 MnGa), 
or even materials with ferrotoroidic ordering (LiCoP0 4 , Co 3 B 7 Oi 3 Br(I)). Related materials 
with additional degrees of freedom include the recently discovered single phase multiferroics 
displaying magnetoelectricity (BiMn0 3 , TbMn0 3 , HoMn0 3 ) where magnetization and po- 
larization are coupled to the lattice and where noncollinear structures arise. The stable 
crystal structure in martensites at high temperatures, often of cubic symmetry, is identified 
as austenite. Near the transformation temperature T c , the coordinated motion of all atoms 
spontaneously breaks the symmetry of the high-temperature phase and several variants of 
the low-temperature phase, the martensite, ensue. This phase is a combination of all the 
individual martensite variants, the fraction of which in the microstructure depends on many 
factors such as the ambient temperature, rate of cooling, or applied stress. The typical sym- 
metries of the martensite include tetragonal (Fe-Pd, In-Tl, Ni-Al(Mn), Fe-Ni-C, BaTi0 3 , 
Ni 2 MnGa), orthorhombic (Au-Cd, U 6 Nb, Ti-Ta(Pd), Cu-Ni-Al(Ti), U), monoclinic (Ni-Ti, 
Ni-Ti- Al(Cu,Fe,Mn,Pd), Cu-Al-Zn(Be), Zr 2 CuCo, Pu), or trigonal/rhombohedral (Au-Cd, 
Ni-Ti, Ti-Ni-Al(Fe)). As recognized early by Landau 1,2 ' 3 , an important condition for dis- 
placive phase transformations to occur is that the symmetry group of the martensite is 
a subgroup of the symmetry group of the austenite. The ratio between the number of the 
symmetry operations of the austenite and the martensite then determines the number of dif- 
ferent martensite variants possible in the microstructure below T c . For example, for a cubic 
to monoclinic phase transformation the initial cubic symmetry is described by 48 symmetry 
elements whereas the final monoclinic symmetry by only 4. Below T c , the high-temperature 
cubic symmetry is thus spontaneously broken into 12 variants of the martensite. Because 
all martensite variants have the same energy, the final microstructure below T c is typically 
composed of regions accommodating these individual variants. When different martensite 
variants are brought together to form an interface, it leads to a strain-matched or strain-free 
twin boundary or domain wall, if properly oriented. Otherwise, there exist transition zones 
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at the domain boundaries that contribute extra compressional and shear energies, and it 
may thus be favorable for the structure to form twin boundary dislocations to lower the free 
energy. Similarly, a habit plane between the parent austenite and product martensite vari- 
ants is an invariant strain plane if properly matched and oriented so that there is a twinned 
microstructure in the martensite with rapidly decaying strain fields in the austenite. Upon 
cooling from the austenitic phase, the individual martensite variants may form a precursor or 
tweed-like microstructure just above the transition temperature which then transforms into 
a fully twinned pattern below the transformation temperature. These morphologies have 
been well characterized experimentally, for example, in the Llo 71-phase of the intermetallic 
Fe-Pd 4 ' 5 ' 6 , Ni-Ti 7 < 8 , Au-Ti 9 , U-Nb 10 , and Cu-Zn-Al 11 . 

Energy minimizing principles are now widely employed to study equilibrium microstruc- 
ture and evolution of martensitic phase transformations. The approach pioneered by Barsch 
and Krumhansl 12 utilizes a nonlinear free energy together with strain inhomogeneity in 
terms of appropriate order parameter strains that drive the transformation from austen- 
ite to martensite. The one-dimensional interface solutions were shown to be soliton-like 
and part of the motivation was to demonstrate how martensite formation can be described 
within a continuum framework without the need for invoking dislocations. These ideas 
were subsequently extended to two dimensions, however, the solutions for the martensitic 
structure were always in terms of displacement fields 1 '' rather than strains. The effects of 
compositional fluctuations and a description of the precursor, the tweed microstructure, in 
Fe-Pd was considered in a Monte Carlo study of Kartha et al. 14 in which the free energy was 
written in terms of strains that contained elastic signatures associated with compositional 
effects. The strains were written as gradients of the displacement field and the Saint- Venant 
compatibility constraint served as the integrability condition for strain fields. Incorporating 
this constraint leads to long-range interactions in the order parameter fields and the mini- 
mum of the free energy is obtained for a twinned microstructure typical of the martensite 
phase. The multiscale consequences of the strain-only model were investigated by Shenoy 
et al. 15 and the nature of the repulsive potential associated with the microstructure, strain 
order parameter dynamics and extensions to two-dimensional ferroelastic transformations 
are reviewed by Lookman et al. 16 Applications to phase transformations in shape memory 
polycrystals and dynamic strain loading in martensites undergoing a cubic to tetragonal 
transformation were further studied by Ahluwalia et al. 1 ' The validity of the Saint- Venant 
compatibility constraint guarantees that the strain field can be obtained from a known dis- 
placement field by taking gradients. On the other hand, if one knows the strain field, the 
displacement field is determined up to a rigid body motion by integration. These statements 
are only true when the material does not contain any topological defects. For example, if 
dislocations are present the displacement field becomes multivalued and its gradient, i.e. the 
strain field, is not defined. Hence, previous studies using a strain-only description are valid 
for defect-free media only. 

Mesoscopic studies of the collective phenomena associated with defects have focused 
primarily on the mechanisms involved in dislocation pattern formation. Groma et al. 18 ' 19 
and Bako et al. 20 formulated a statistical model for the evolution of the dislocation den- 
sity in isotropic bodies, where the stress field associated with each dislocation is given 
analytically" 1 '"^. For the case of single slip in two dimensions, the dislocation density evolves 
according to a Fokker-Planck equation in which the Peach-Koehler force" 5 on each disloca- 
tion is determined from the known stress field. A different framework, based on statistical 
studies of dislocation patterning developed by El-Azab 24 ' 25 , makes closer connection with 



3 



the Kroner's 26 continuum theory of dislocations. It accounts not only for long-range interac- 
tions between dislocations but also treats each discrete slip system separately. By calculating 
pair correlations between dislocations in this statistical model, Zaiser et al. 2 " demonstrated 
that dislocation systems exhibit a patterning instability and this leads to the formation of 
dislocation walls perpendicular to the glide plane. These statistical ensembles of disloca- 
tions have been shown to exhibit intrinsic spatio-temporal fluctuations with scale-invariant 
characteristics, long-range correlations and emergence of strain bursts (for a recent review, 
see Zaiser 28 ). In addition, continuum theories of dislocations and self-stresses as developed 
by Kroner 26 , Kosevich 2!} and others have been applied to studies of dislocation patterning. 
Among the most prominent are contributions of Kratochvfl and Sedlacek 30 and Sedlacek 
et al. 31 where the evolution of the dislocation density is formulated in both the Eulerian and 
the Lagrangian frames. A phase field formulation of the dislocation patterning in isotropic 
media was developed by Rickman and Vinals 32 , where the dislocation density tensor is 
obtained by minimizing the free energy. However, this model does not include coupling be- 
tween the dislocation density and the underlying crystal structure and, therefore, it does not 
exhibit any structural phase transition. Most recently, Roy and Acharya 33 implemented the 
Kroner's theory to study dislocation patterning using the Finite Element Method (FEM). 
A similar approach was adopted by Limkumnerd and Sethna ,l to formulate a mesoscopic 
Landau theory in which the free energy is written in terms of the plastic distortion tensor. 
The plastic distortion field that minimizes the free energy is then used to calculate the dis- 
location density which plays the role of the order parameter. The added advantage of this 
formulation is that the same model applies to dislocation motion by pure glide and by a 
combination of glide and climb and, therefore, it allows for studies of dislocation patterning 
in isotropic materials at both low and high temperatures. However, the stress associated 
with individual dislocations is still calculated using the isotropic elasticity 21 ' 22 . 

Our objective in this paper is to incorporate dislocations into the Landau theory to study 
martensitic phase transformations in materials containing defects. Unlike the previous stud- 
ies cited above, we consider an anisotropic medium that is described by the elastic constants 
corresponding to the high-temperature cubic phase. Utilizing Kroner's 2 '' continuum theory 
of dislocations, we show that the presence of dislocations induces incompatibility between 
the elastic components of the strain tensor field, and this is connected to the gradients of 
the dislocation density. The presence of dislocations is responsible for a nonlocal coupling 
of the incompatibility field with the order parameter and, as a consequence, the evolving 
martensitic texture is affected by the finite density of dislocations. Minimizing the free 
energy subject to the incompatibility constraint for a given distribution of dislocations gen- 
erates a stress field that corresponds exactly to this distribution of defects. By inserting a 
single edge dislocation into an otherwise ideal crystal we show that the order parameter field 
that minimizes the free energy subject to this incompatibility constraint yields the correct 
long-range stress field around this dislocation. The fact that the stress field in a generally 
anisotropic material with arbitrary distribution of dislocations can be calculated by merely 
minimizing the free energy means that we can easily calculate the Peach-Koehler forces that 
act on these individual dislocations. These forces are then used in the Fokker-Planck equa- 
tion for an evolution of the dislocation density. The procedure outlined above represents a 
self-consistent scheme that is solved recursively. In the first step, the order parameter field 
is calculated by minimizing the free energy subject to the incompatibility constraint that 
is obtained from the given dislocation density. For the known order parameter field, the 
Peach-Koehler forces on individual dislocations are calculated and utilized in the evolution 
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equation to update the dislocation density and thus also the corresponding incompatibility 
field. 

The mesoscopic nature of our approach makes this work different from the more 
nanoscale-based models, in particular, the phase field microelasticity theory of Khachatu- 
ryan 35 and Wang et al. 36 . In their work the dislocation loops are viewed as coherent platelet 
inclusions that expand, interact with other loops in the same and other slip systems and 
annihilate in response to their internal long-range strain fields and externally applied load. 
The dislocation content of each slip system is described by an integer- valued density function 
rj that specifies the number of dislocations with prescribed direction of the Burgers vector. 
The total energy is written in terms of stress-free strain due to individual dislocation loops 
and the evolution of the order parameter field r\ is studied using a Langevin dynamics based 
on the time-dependent Ginzburg-Landau equation. This model has been successfully uti- 
lized in studies of heterogeneous nucleation of martensite in the parent austenitic martix 
(for an excellent review, see Malygin 3 ' ). The martensitic embryo was shown to grow in- 
side the dislocation loop that expands in response to the external load, giving rise to the 
so-called stress-induced martensite. In contrast, within the Barsch-Krumhansl formulation 
of martensite 12 , which is utilized in this paper, twins nucleate spontaneously by lowering 
the temperature and without the need of an externally applied load provided the system 
contains a certain minimum degree of strain inhomogeneities induced by thermal fluctua- 
tions of the lattice. In the present work, this is no longer needed because the heterogeneous 
nucleation of the martensite takes place readily on preexisting dislocations. 

The plan of our paper is as follows. In Section II we review the continuum theory of 
dislocations and show how we include the dislocation density as the source of the incompat- 
ibility. In Section III we eliminate this constraint in favor of long-range interactions in the 
dislocation density and its nonlocal coupling with the elastic strains. In Sections IV and V 
we show how the total free energy is minimized using relaxational dynamics and how the 
dislocation density is evolved using the Fokker-Planck equation, respectively. In Section VI 
we utilize the theory developed in this paper to study martensitic phase transformations 
and dislocation patterning in single crystals of Fe-30at.%Pd alloys. In the first case study 
we consider that the dislocation density is fixed and thus the free energy is minimized solely 
by the order parameter field. We demonstrate that the increase of the dislocation density 
induces long-range internal strains in the material, and these give rise to stress-induced 
martensite even above the temperature T c for a defect-free material. Below T c the morphol- 
ogy of the martensite changes from well-defined (110) twins at low dislocation density to a 
twin-free order parameter field at high dislocation density. In the second case we consider 
ideally mobile crystal dislocations. We show that below T c the order parameter field evolves 
into a series of twins corresponding to the two variants of the martensite. The equilibrium 
dislocation density is characterized by regions of high dislocation density in which the dislo- 
cations are arranged into walls at twin boundaries, separated by regions of low dislocation 
density. We show that the neighboring dislocation walls are formed by opposite Burgers 
vectors perpendicular to the twin boundaries. These correlated domains have been observed 
not only in molecular dynamics and Monte Carlo simulations 19 ' 38 , but also in experiments 
on Ag''' } , Ni-Ti 7 and Fe-Pd 40 . 
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II. CONTINUUM THEORY OF DISLOCATIONS 



In order to demonstrate how the presence of defects breaks the single-valuedness of the 
displacement field, we will first perform the following thought experiment as proposed orig- 
inally by Kroner 26 ' 41 . Consider a macroscopic single crystal whose structure belongs to a 
well-defined space group and its unit cell is defined by the lattice parameters a = {a , b , c } 
and angles ao = {ao,/?o ; 7o}- Only those microscopic details associated with the crystal 
structure that manifest themselves at the mesoscopic level will be taken into account. This 
leads us to discretize the medium into a finite number of mesoscopic cells with lattice param- 
eters a = {a,b,c} and angles a = {a, (3, 7}, where each such cell includes a finite number 
of the crystallographic unit cells. The dislocation content of each mesoscopic cell can thus 
be characterized by the so-called net Burgers vector B that is defined as a vector sum of 
the Burgers vectors 6; of the crystal dislocations in all embedded crystallographic unit cells. 
This is shown schematically in Fig. 1. Each of these "net dislocations" characterized by 
its Burgers vector B, causes a certain plastic distortion (3 P that is proportional to the mag- 
nitude of B. We require that the continuity of the body be maintained for any distribution 
of the net Burgers vectors, i.e. for any corresponding plastic distortion (3 P . In other words, 
by inserting dislocations in the originally defect-free medium we are not allowed to create 
microcracks that would destroy the continuity of the body. This requirement is equivalent 
to demanding that the total distortion field (3 l be curl-free, i.e. 



V x (3* 







(1) 



Here, (3 = (3 + (3 P is the total plastic distortion 50 that is written as a sum of its elastic 
part (3 and the plastic part (3 P . Each of these distortions can be expressed as a sum of their 
symmetric parts, i.e. strains £*, e, e p , and antisymmetric parts corresponding to rotations 
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FIG. 1: Each mesoscopic cell (solid square) is comprised of a finite number of crystallographic unit 
cells (dotted gray squares). The net Burgers vector B is a vector sum of the Burgers vectors of 
individual crystallographic dislocations (gray arrows). 

In general, two kinds of plastic distortions of a mesoscopic cell can take place depending 
on the way the adjacent cells are distorted and these are shown in Fig. 2. If the net Burgers 
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vector in a given cell is the same as those in the adjacent cells, all these cells are distorted the 
same way and, therefore, the continuity of the body is locally preserved. Hence, V x (3 P = 
and, with the help of (1) and /3* = (3 + (3 P , this means that the elastic part of the strain 
field is also curl-free, i.e. V x (3 = 0. In this case, the above-mentioned requirement of the 
continuous medium leads to the well-known Saint- Venant elastic compatibility constraint ' 1 



V x V x e = 



(2) 



Since the elastic strains are compatible in the sense of the Saint- Venant law, this plastic 
distortion is referred to as compatible"' 2 . In contrast, one can imagine a more general case 
where the adjacent cells are characterized by different net Burgers vectors, which means 
that the cells are distorted differently. Hence, the plastic part of the distortion tensor is no 
longer curl- free and, instead, — V x f3 p = a, where a is the tensor of the density of net 
Burgers vectors 53 . If this plastic distortion acted alone, it would cause disregistry between 
neighboring mesoscopic cells 2 ' 1 and thus contradiction of the requirement that we set forth by 
(1). However, cohesive forces of the matter act to remove this disregistry and this relaxation 
proceeds purely elastically. Clearly, in order to satisfy (1), the elastic strain field cannot 
be arbitrary but has to satisfy the constraint V x (3 = a. Performing the curl of this 
equation and taking its symmetric part then leads to an incompatibility constraint between 
the components of the elastic strain tensor, 



V x V x e = r) , 
where r) is the so-called incompatibility tensor defined as 

77 = sym (Vxq) 



(3) 



(4) 



In this case, the individual components of the elastic part of the strain tensor, e, are not 
compatible in the sense of the Saint- Venant law (2). Hence, this plastic distortion is called 
incompatible and the degree of this incompatibility is quantified by the symmetric tensor 77. 
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(a) compatible plastic distortion 
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(b) incompatible plastic distortion 



FIG. 2: Two kinds of plastic distortion of two neighboring mesoscopic cells. The compatible plastic 
distortion (a) maintains the continuity of the body, while the incompatible plastic distortion (b) 
requires elastic relaxation to keep the body crack-free. 



It is important to emphasize that since each net Burgers vector is a sum of many Burg- 
ers vectors of the underlying crystal dislocations, this coarse-graining procedure essentially 
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determines the vector corresponding to the excess of the crystallographic Burgers vectors. 
Moreover, since each mesoscopic cell contains a large number of crystal dislocations, the 
density of net Burgers vectors can be approximated by a tensor field that is continuous 
throughout the entire space. In the Cartesian coordinate system with axes x±, x 2 , x$, this 
density is represented by the Nye 12 tensor a^, the components of which relate to the net 
dislocations with line directions parallel to the Xj axis and the Burgers vectors parallel to the 
Xj axis. Hence, the diagonal elements of at correspond to screw components, while the off- 
diagonal elements to edge components of the net Burgers vectors. One can thus determine 
the density of the net Burgers vectors as ay = Bj/Si, where Bj is the j-th component of 
the net Burgers vector and Si are the components of the vector normal to the oriented area 
of the mesoscopic cell pinned by the dislocation line. The tensor cn^ should not be confused 
with the dislocation density p that is, by definition, the total length of all dislocation lines 
that populate the medium divided by its volume. 

For simplicity, consider now a two-dimensional plane strain problem in which the only 
nonzero components of the elastic strain tensor are en, E\ 2 , and £22- Hence, it is straight- 
forward to prove that the only component equation of (3) that is not identically zero is 

fen - 2d 12 e 12 + <9ii£ 2 2 = V33 , (5) 

where d\j = d 2 /dxidxj. It is important to recognize that the scalar incompatibility field 
7733 is nonzero wherever the distribution of dislocations in the medium causes the plastic 
distortion to be incompatible in the sense explained above. Hence, the relation (5) represents 
a constraint that the elastic strains have to satisfy in order to maintain the continuity of 
the body that has been broken locally by the incompatible plastic distortion. From (4), the 
incompatibility field r/33 can be determined as 

7733 = dxa Z2 - <9 2 «3i , (6) 

where di = d/dx-i. Since only the components 031 and a^ 2 appear in this expression, 
only edge dislocations with their line directions parallel to £3 and the Burgers vector com- 
ponents along the x\ and x 2 axes contribute to the elastic strain incompatibility in this 
two-dimensional case. 

In order to develop a clear link between the microscopic crystal dislocations and their 
coarse-grained mesoscopic manifestation by the net Burgers vector B, we will now consider 
a finite number of discrete slip systems s. In each such system, the Burgers vector of crystal 
dislocations b s is known, e.g. 6 s = 1/2(110) (in units of the lattice parameter do) for edge 
dislocations in face-centered cubic crystals. To each mesoscopic cell we can then attribute 
the net Burgers vector 

B(r) = J2N s (r)b s , (7) 

s 

where N s is the number of crystal dislocations in the slip system s with the Burgers vector 
b s . Instead of working with integral values N s , it is convenient to define a number density 
of crystal dislocations in the slip system s as n s = N s /S ce u, where S ce u is the area of a 
mesoscopic cell, and regard this as a continuous variable. Hence, the density of net Burgers 
vectors can be written as 

a 3i (r)=B l (r)/S cell = Y,n s (r)b s i • (8) 
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Substituting this expression in (6) yields the strain incompatibility in terms of the density 
of crystal dislocations in individual slip systems: 

where is the Levi-Civita tensor. The expression (9) represents a "recipe" for coarse- 
graining the density of crystal dislocations n s in individual slip systems s into a continuously 
distributed incompatibility field 7733. It is important to emphasize that each slip system s 
contains crystal dislocations with positive and negative Burgers vectors of the same magni- 
tude. The dislocation density n s can thus be written as n s = n s+ — n s ~, where n s+ and n s ~ 
are non-negative densities of crystal dislocations with positive (b s ) and negative (—b s ) Burg- 
ers vectors, respectively. In Section V we will show that this distinction between positive 
and negative Burgers vectors of crystal dislocations is required for a systematic evolution of 
the dislocation density. 



III. MESOSCOPIC FREE ENERGY 

Consider an elastically anisotropic body of which every element is subjected to a generally 
nonuniform stress tensor cr and the response to this particular loading is characterized in 
each such element by the elastic strain tensor e. If we consider that a linear relation between 
the applied stress and induced strain applies, the free energy of this deformed medium equals 
the total strain energy, i.e. F = J v ^(JijEijdr. Writing = c^kiSku where c^ki is the elastic 
stiffness tensor, one arrives at the free energy 43 

F = J ^CijkiEijEki&r . (10) 

The elastic stiffness tensor has generally 21 independent elastic constants but any sym- 
metry of the underlying crystal structure reduces this number. If the stress tensor is 
written in the Voigt notation as <X; = [cr n <7 2 2 033 °23 013 o- 12 ] T and the strain tensor 
as Sj = [en £22 £33 £23 £13 Zu} 7 ', the free energy (10) can be expressed equivalently as 
F = J v ^CijEiEjdr, where is the (6 x 6) symmetric elastic stiffness matrix. For simplic- 
ity, we will be concerned in the following with cubic symmetry and this is characterized by 
three independent elastic constants Cn, C12, and C 44 . The corresponding free energy for 
cubic symmetry then reads 

F = J v ^Cu(el 1 +E 2 22 ) + C 12 E 11 e22 + 2Cuel 2 ^dr , (11) 

where we write the strains again in their usual two-index notation. 

For simplicity, we will specialize in the following to the case of the square to rectangle 
phase transformations 51 that can be thought of as a two-dimensional reduction of the cu- 
bic to tetragonal (or tetragonal to orthorhombic) phase transformation frequently observed 
in shape memory alloys such as Fe-Pd, In-Tl, Ni-Al(Mn), ternary alloys Fe-Ni-C, magne- 
toelastic alloy N^MnGa, or even perovskites such as BaTiC>3 that exhibit strain-induced 
polarization. A naive way to identify a phase transformation would be to define the order 
parameter as a ratio of two lattice parameters. However, this ratio alone would not dis- 
tinguish between cubic and rhombohedral symmetry and, therefore, one order parameter is 
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generally not sufficient to identify the crystal structure. In the following, we consider three 
scalar order parameter fields 11 , 



ei 



V2 



£ll + ^22 J 



e 2 



V2 



En — £22j 



e 3 = £12 



12) 



where ei measures the isotropic dilation, e2 the deviatoric change of shape, and e 3 the change 
of the right angle caused by the shear. These three fundamental modes of deformation of 
an element in the body are shown in Fig. 3, where €2 = corresponds to the austenite and 
negative/positive to the two variants of the martensite. Here, 62 serves as the primary 
order parameter for the square to rectangle transformation, whereas e\ and e 3 are secondary 
order parameters. In terms of these fields, the Hookean elastic free energy (11) for this 
two-dimensional problem is 



A x 



e? + 4^ + 



.4 



3„2 



dr 



(13) 



where the coefficients are related to the elastic constants 55 as Ai 
and A 3 = 4C44. 



C11 + C12, A 2 — C\\ — C 



12, 





ei 



e 2 



f'3 



FIG. 3: The three order parameters used to study the square to rectangle phase transformation in 
crystals of cubic symmetry; red=positive, blue=negative value. 

For studies of first order phase transitions the order parameter €2 is expanded in even 
powers up to the sixth order as odd powers are not allowed by symmetry. In addition, it 
is customary 14 ' 15 ' 16 ' 44 to incorporate a gradient term proportional to |Ve2| 2 that represents 
the energy cost for spatial variation of the order parameter, and the strain energy due to 
coupling of the internal strain with the externally applied stress field. Consequently, the 
free energy can be written as 



F = / {jfoc[e 2 ] + fnonioc[ei, e 3 ] + / gr ad[e 2 ] - fioad[ei, e 2 , e 3 ]} dr 
Js 

where the various energy densities are 

f nordoc\&\i 63] 
fgrad[^2\ 

fioad[ei, e 2 , e 3 ] 



(14) 



.<4 2 2 ^4 ^6 



A x 



A 



3 „2 
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|Ve 2 | 



(15) 
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It is important to emphasize that the strain energy density as given by fi oa d applies only to 
the linear-elastic case, i.e. when the internal strains induced in the body by the external 
loading vary linearly with this applied stress. 

In defect-free media the individual components of the elastic strain tensor e are related by 
the Saint- Venant constraint (2) and this guarantees that the strains can be determined from 
the known displacement field by taking its gradients. In contrast, the presence of dislocations 
causes discontinuities in the displacement field and, therefore, the integrability condition no 
longer applies. Hence, the strains are incompatible in the sense of the Saint- Venant condition 
and this is expressed by (3), where 77 represents the "strength" of this incompatibility. In the 
case of plane strain in the (xi,x 2 ) plane the incompatibility constraint (5) can be expressed 
using (12) in terms of the order parameters as 

V 2 ei - (<9n - <9 22 )e 2 - VSd 12 e 3 = r) 33 \fl . (16) 

Now, suppose that we know the density of crystal dislocations n s in each slip system s. 
Hence, the corresponding incompatibility field 7733, obtained from (9), is known at every 
point. We then seek the fields e±, e 2 , e 3 that minimize the free energy (14) subject to the 
incompatibility constraint (16). The additive nature of the free energy allows us to perform 
this minimization in steps. Firstly, we calculate the fields e\ and e 3 that are constrained 
by the incompatibility condition (16) and, in the second step, we obtain the primary order 
parameter e 2 by minimizing the free energy. The part of the free energy that depends 
explicitly on e x and e 3 reads 

^13= / { /nonioc [ei, e 3 ] + /w[ei,e 2 ,e 3 ] + AG [ei,e 2 ,e 3 , 7733]} dr , (17) 
Js 

where A is the Lagrange multiplier that incorporates the incompatibility constraint G = 
V 2 ei — (<9n — <9 22 )e 2 — V8d\ 2 e 3 — r] 33 ^/2 = 0. The fields e\ and e 3 that minimize Fi 3 are then 
obtained from the stationary conditions 56 5Fi 3 /5ei = 0, 5F 13 /5e 3 = 0, and 5Fi 3 /5\ = 0. 
When evaluated in k-space, these conditions provide closed-form expressions for the sec- 
ondary order parameter fields ei, e 3 in terms of the incompatibility 7733 and the components 
of the externally applied stress tensor a if 

ei (k) = Q t (k)e 2 (k) - V2R t (k) V33 (k) - SMa^k) - T t (k)a 12 (k) , (18) 

where i — 1,3. Here, Q x ~ — ky, Q 3 ~ k x k y {k 2 x — k 2 y ), Ri ~ k 2 , R 3 ~ k x k y) S\ ~ k 4 , 
S 3 ~ k x k y k 2 , Ti ~ k x k y k 2 , and T 3 ~ k 2 k 2 are k-space kernels that we write explicitly in 
Appendix B, and ajj = an + cr 22 . Eq. (18) represents the most general form " valid for 
the plane strain case in which the stress tensor has all components nonzero. With the help 
of (18), the nonlocal part of the free energy in (15) can be written as a functional of the 
primary order parameter e 2 , the incompatibility 7733, and the components of the stress tensor 
(Jij only: 

fnomoc(k) = M^l[e 2 (fc)] 2 - V2B 13 (k)e 2 (k) V33 (k) - S j4Q (fc)e 2 (fc) + 

+C u (k)[ V33 (k)] 2 + V2^ AR (k) V33 (k) + ^^ , (19) 

where the k-space kernels S^q, T>ar and E^, written explicitly in Appendix B, depend on 
the external stress field. In (19), the first term represents nonlocal interactions in the e 2 



11 



field, whereas the second and third terms are couplings of e 2 with the incompatibility 7733 and 
the external stress <Jij, respectively. The fourth and the fifth terms are contributions from 
the incompatibility and its coupling to the external stress field, and the last term represents 
the shift of the free energy by the external stress field. If the incompatibility 7733 vanishes 
and no external stress is applied, only the first term remains and we recover the nonlocal 
expression for an unloaded defect-free medium 14 ' 16 . 

In a similar way as we expressed f non i oc in terms of e 2 , 7733, and only, we can utilize 
(18) to obtain a reduced expression for the strain energy density fi oa d- For plane strain, a 
completely general stress state leads to fi oa d = \o\\E\\ + 0^2^12 + \^n^22- Expressing the 
strains in terms of the order parameters ei, e 2 , e$, and using (18), one arrives at the strain 
energy density 

fioad(k) = W Q (k)e 2 (k) - V2W R (k) V33 (k) - Wx(k) , (20) 

where the kernels Wq, Wr and Wt, are again written explicitly in Appendix B. The remain- 
ing two free energy densities, i.e. f\ oc and f gra d, are functionals of e 2 only and are determined 
uniquely by (15). 



IV. RELAXATION OF THE PRIMARY ORDER PARAMETER FIELD 

Since all constituents of the free energy are now functionals of the primary order parame- 
ter e 2 , the incompatibility 7733 and the stress tensor cr^ and, assuming that the incompatibility 
and the external stress field change slowly relative to e 2 , i.e. they remain approximately 
constant on the time scale of relaxation of e 2 , it is straightforward to find the field e 2 that 
minimizes the free energy. This minimization, i.e. the solution of the equations SF/Se2 = 
and 8 2 F/8e\ > 0, cannot be performed analytically. However, we may formulate a relax- 
ational dynamics 58 for e 2 that will follow the path of the steepest descent of the free energy 

where T plays a role of the mobility parameter. Writing (21) as a difference scheme, i.e. 
e 2 (t + At) = e 2 (t) — AtT 5F/5e 2 , it is clear that T merely renormalizes the time step At. 
Moreover, since 5F/5e 2 = d{f ioc + fnonioc + fgrad - fioad) /de 2 , the right-hand side of (21) can 
be calculated easily by taking derivatives of the previously derived free energy densities: 



dfiocjr) 
de 2 {r) 



A 2 e 2 (r) + J B[e 2 (r)] 3 + C[e 2 (r)] 5 



= ^3(fc)e 2 (fc) - V2B 13 (k) V33 (k) - S AQ (fc) , 

= _^ 2 V 2 e 2 (r) , (22) 

<9e 2 (r) 

dfload(k) 

de 2 (k) 



W Q (k) . 



If one considers a defect-free medium and no external stress is applied, the second and the 
third terms on the right-hand side of df non i oc / de 2 vanish and also dfi oa d/ de 2 = 0. In this case, 
the nonlocal free energy density (19) reduces to f non ioc = (Ai3(k)/2) [e 2 (fc)] 2 which is identical 
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to the form obtained by Kartha et al. 14 Recognizing that the kernel A\^(k) ~ (k 2 — k 2 ) 2 is 
minimized when k x = ±k y , it directly follows that the system minimizes its free energy by 
aligning the nonzero components of the order parameter field €2 along any of the two k-space 
diagonals, as shown in Fig. 4a. The system thus develops diagonal striations in real space, 
similar to the tweed microstructure that is a general feature of many martensites 1 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 . 




(a) 



(b) 



FIG. 4: Schematic illustration of the texture in ei that minimizes the first two terms in the 
nonlocal part of the free energy density (19). (a) minimizes Ai^{k)[e2{k)] 2 - the two orientations 
of twins yield the same free energy; (b) minimizes — Bi^{k)e2{k)r]^{k) - gray squares correspond 
to e2(k)r]3s(k) > and black squares to e2(fc) 7 ?33 (fc) < 0. 



The presence of dislocations changes the evolving order parameter texture considerably 
and this can be demonstrated by examining the first three terms in (19). As mentioned 
above, the first term is minimized when all nonzero components of C2{k) are aligned along 
the k-space diagonal, as shown in Fig. 4a. In the second term, the kernel — 5 13 (fc) ~ —k 2 +k 2 
is minimized when k = (±7r//,0) and maximized when k = (0, where I is the real- 

space width of the simulated domain. Therefore, the second term in (19) is minimized when 
positive e2(k)rj 33 (k) is at the minimum of — B\s(k), i.e. at k = (±7r//, 0) and, simultaneously, 
negative e2(fc) ? ?33 (k) is at its maximum, i.e. at k = (0, ±tc/1). This is shown in Fig. 4b. 
In the real space the minimum of the second term in (19) thus corresponds to a three-state 
"checkerboard" pattern with the periodicities along the x± and £2 axes equal to twice the 
width of the mesoscopic cell. The size of the mesoscopic cell thus imposes an intrinsic length 
scale in the order parameter field. Similarly, one can identify the pattern that minimizes the 
third term in (19) which now depends on the external stress field. Since the first three terms 
in fnonioc cannot be minimized simultaneously when a finite incompatibility is introduced, 
they naturally compete with each other. Hence, the minimization of the free energy by the 
field €2 subject to a fixed distribution of dislocations is inherently frustrated and €2 does not 
always evolve into a well-defined diagonal texture as it does in defect-free materials. We will 
see later that the diagonal texture is preferred at low dislocation densities, whereas at high 
dislocation densities the terms containing the incompatibility 7733 become significant and the 
texture in e 2 tends to that corresponding to Fig. 4b. 



V. EVOLUTION OF THE DISLOCATION DENSITY 

From the order parameter field €2 that minimizes the free energy, we can obtain ei and 
using (18) and subsequently en, E\2 and e 2 2 using (12). Assuming linear-elastic dependence 
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between the stresses and strains and using the three anisotropic elastic constants Cn, C\ 2 
and C44, the components of the internal stress field are obtained from: 
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^33 (r) = ~ [o-n(r) + <T 22 (r)] . 

Note that just by minimizing the free energy subject to the elastic strain incompatibility 
constraint we obtained an internal stress field that is a superposition of the elastic stress fields 
of individual net dislocations. In order to demonstrate this, we show in Fig. 5 the calculated 
stress field around one edge dislocation in an isotropic medium with its Burgers vector 
along the positive x\ direction. The overall distribution of stresses is in excellent agreement 
with the formulas derived from isotropic elasticity (see Hirth and Lothe 21 ). However, in 
the derivation of the free energy no constraints were imposed on the anisotropy of elastic 
constants and, therefore, the stress fields of the dislocations in arbitrarily anisotropic media 
can also be calculated just by minimizing the free energy. 




(a) <7n(r) (b) cr 22 (r) (c) a 12 (r) 



FIG. 5: Stress field around one edge dislocation with its Burgers vector parallel to the x\ axis 
(blue=negative, red=positive values). Because periodic boundary conditions are used in this 
calculation, the stress field shown here corresponds to one dislocation from a periodic array of 
dislocations of the same kind. 

It is important to realize that each slip system contains crystal dislocations with positive 
and negative Burgers vectors and in our mesoscopic description we have to treat the densities 
of these dislocations separately. The explanation is provided in Fig. 6 where we show 
schematically two different sizes of the mesoscopic cell and the corresponding parallel slip 
planes (dotted lines) that pass through these cells. If the mesoscopic cell is of the same size 
as the crystallographic unit cell (Fig. 6a), only one slip plane corresponding to each slip 
system s passes through this cell. In this case, positive and negative dislocations meeting in 
this cell annihilate and the only relevant quantity is the excess density n s . However, if the 
mesoscopic cell comprises a number of unit cells, shown schematically in Fig. 6b, which is also 
the case in our model, several parallel slip planes pass through the cell. If we consider that the 
dislocations move by pure glide, the mesoscopic cell can thus contain simultaneously positive 
and negative dislocations in parallel slip planes and only a limited number of dislocations can 
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annihilate. Clearly, if we are to reproduce correctly the dislocation content of the mesoscopic 
cell we have to consider the densities n s+ and n s ~ separately. 
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FIG. 6: Schematic illustration of the slip planes (dotted lines) corresponding to a particular slip 
system s passing through the mesoscopic cell (black square) when it comprises only one (a) and 
many (b) crystallographic unit cells (gray squares). 



Without the presence of dislocation sources, the densities n s+ and n s ~ of crystal dislo- 
cations in each slip system s integrated through the simulated domain must be conserved 
individually. This implies that the evolution equation for these dislocations densities must 
take the form of a continuity equation. Since we consider here that crystal dislocations move 
only by glide in their corresponding slip planes, the evolution equation for the dislocation 
densities reads 

^-^•[^(r.rtt)] . (24) 

Here, F s g f ide (r, t) represents the glide component of the Peach-Koehler force on the crystal 
dislocations with the densities n s± in the mesoscopic cell at r and time t, and can be 
calculated as follows. If the internal stress tensor is known at r, the components of the total 
Peach-Koehler force on each crystal dislocation with the Burgers vector ±6 S within the same 
mesoscopic cell can be calculated 21 ' 22 as = ^fejk&jibf ■ The glide force is then determined 
by projecting this force into the corresponding slip plane, i.e. F s lide = (F s± ■ e s± )e s± , where 
e s± = b s /|6 s |. The dislocation densities can then be updated using (24). For each slip 
system s the fields n s+ and n s ~ then enter (9) to calculate the incompatibility field 7733. 
One can thus recalculate the order parameter that minimizes the free energy (14) subject 
to this updated distribution of incompatibilities. The calculation is regarded as complete 
when the free energy is minimized by the field e% and, simultaneously, the corresponding 
stress field does not induce significant changes in the dislocation densities n s+ and n s ~ . 
This process thus represents a self-consistent procedure for the simultaneous evolution of 
the primary order parameter field €2 and of the densities of crystal dislocations in individual 
slip systems. 



VI. SIMULATIONS 



The material considered in these simulations is the shape memory alloy single crystal 
of Fe-30 at.%Pd for which the temperature dependence of the elastic constants, measured 
by Muto et al. 45 , is parametrized together with other coefficients entering the free energy 
functional in Appendix C. In our calculations the orientation of the crystal is chosen such 
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that the x± axis coincides with the [100] direction, and the X2 axis is parallel to the [010] di- 
rection. Each mesoscopic cell is characterized by 1000 x 1000 crystallographic unit cells with 
the lattice parameter ao = 3.8 A. Hence the width of each mesoscopic cell is a = 0.38 pm 
and its area S ce u ~ 0.14 /1m 2 . The initial values of the order parameter e 2 in each meso- 
scopic cell are drawn randomly from a uniform distribution with limits (—0.04, 0.04) which 
include the values of e 2 that correspond to the minima of fi oc below T c . The simulated two- 
dimensional domain contains 128 x 128 mesoscopic cells with periodic boundary conditions 
that are automatically imposed by the k-space calculations. The width of the simulated 
domain is 48.64 pm. For simplicity, no external stress is applied in any simulation. The 
minimization of the free energy is regarded as complete when bF jbe-i in (21) becomes less 
than 0.001 eV/A 3 . 

The TEM observations of dislocations in polytwinned Fe-Pd thin foils due to Xu et al. 6 
reveal that the active mode accommodating plastic strain in this material is due to the glide 
of 1/2(110) dislocations. Hence, in our two-dimensional simulations we will consider only two 
slip systems s, populated by crystal dislocations with Bugers vectors 6 1± = ±1/2[110] and 
b 2± = ±1/2 [1 10]. Each mesoscopic cell is initially assigned finite densities of these crystal 
dislocations, n s+ and n s ~ , each of which is chosen at random from a uniform distribution 
between and n max , where n max is varied to arrive at different dislocation densities. For 
each simulation the actual dislocation density p is given below. 

A. Martensitic texture as a function of the density of dislocations 

We first study the distribution of the order parameter e 2 that minimizes the free energy 
(14): (i) in dislocation-free material, and (ii) subjected to fixed dislocation density. Three fi- 
nite densities of crystal dislocations are considered, given by p = {6x 10 14 , 10 15 , 8x 10 15 } m -2 . 
The free energy (14) is then minimized subject to the distribution of incompatibilities 7733 
that are derived from the given dislocation density using (9). No evolution of the dislocation 
density is allowed in this case which corresponds to a hypothetical situation where all dislo- 
cations are immobile. This minimization is performed for the temperatures both above T c , 
where the austenite is stable in the defect-free medium, and below T c , where the martensite 
is the stable phase. 

In Fig. 7 we show how the order parameter field €2 that minimizes the free energy (14) 
chang GS £lS db function of the density of dislocations. The upper row corresponds to the 
temperature above T c whereas the lower row to the temperature below T c . Above T c , the 
system is progressively driven away from its free energy minimum (e 2 = 0) the higher 
the dislocation density p. This results in stabilization of the martensitic phase above T c 
in the order parameter field e 2 without a well-defined texture (Fig. 7b-d). Below T c , the 
twinned microstructure that minimizes the free energy in dislocation-free materials changes 
at finite densities of dislocations as a consequence of the competition of the first two terms 
on the right-hand side of (19). For low densities p, the first term dominates and the order 
parameter field e 2 is characterized by a twinned microstructure. With increasing dislocation 
density, the second term in (19) becomes of the same order as the first and this competition 
gradually causes elimination of the twins (Fig. 7b-d). For large dislocation densities, i.e. 
Fig. 7d, the strain incompatibility 7733 completely dominates the minimization of the free 
energy and the same martensitic texture is obtained both above and below T c for the defect- 
free medium. These results are consistent with our previous observations that dislocations 
alter the character of the martensitic texture (see Fig. 4). 
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T>T C 




(a) (b) (c) (d) 

FIG. 7: Distribution of the order parameter field e2(V) above and below T c for dislocation-free 
medium as a function of the density of dislocations, (a) no dislocations, (b) p = 6 x 10 14 m~ 2 , (c) 
p = 10 15 m~ 2 , (d) p = 8 x 10 15 m" 2 . 

Recall that the results shown in Fig. 7 correspond to a highly idealized case where all 
dislocations are immobile and the dislocation density thus cannot evolve in response to 
changes in the order parameter field. In the following example, we remove this constraint 
to simultaneously evolve the order parameter field and the dislocation density to provide a 
clearer picture of dislocation pattern formation below T c . 

B. Formation of dislocation walls at twin boundaries 

The initial distribution of the order parameter field above and below T c was obtained 
by minimizing the free energy for a dislocation-free material. In each mesoscopic cell, the 
two directions of the Burgers vectors in both slip systems s were then assigned randomly 
a dislocation density between and n max = 10 14 m -2 which yields an actual dislocation 
density p = 2 x 10 14 m -2 . We then calculated the order parameter field that corresponds to 
the given initial distribution of dislocations. This e% field is used to update the dislocation 
density using (24) where all dislocations have the same mobility, i.e. D is the same for 
the evolution of all densities n s± . Since no external stress is applied the evolution of the 
dislocation density is driven entirely by the evolving martensitic texture and the long-range 
interactions between dislocations. 

Above T c the order parameter field e2 is close to zero everywhere and the corresponding 
dislocation density is spatially uniform. Below T c the field e2 is represented by a series of 
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(a) 



(b) 



FIG. 8: The order parameter field e 2 (a) and the distribution of the magnitudes of net Burgers 
vectors \B\ (b) that minimize the free energy (14). The blue and red regions in (a) are the two 
variants of the martensite while the twin boundaries correspond to the austenite phase. In (b) the 
bright spots correspond to regions of high dislocation density and dark spots to low dislocation 
density. The inset of (b) shows the orientations of net Burgers vectors in the region marked by the 
gray square. 

twins corresponding to the two variants of martensite (see Fig. 8a). As can be seen from 
Fig. 8b the scalar density of net Burgers vectors, |-B(r)| = a/ Bi{r) 2 + B 2 (r) 2 , corresponding 
to the minimum of the free energy, is characterized by twin boundary dislocations with the 
Burgers vectors parallel to the [110] and [110] directions. Hence, the twin boundaries are 
decorated by dislocation walls with the net Burgers vectors parallel to the (110) direction, 
as shown in the inset of Fig. 8b. This results in the formation of alternating positive and 
negative dislocation walls along the twin boundaries. Since all (110) twins have the same 
energy, the free energy can also be minimized by forming the 45° (rather than 135°) texture 
in e 2 (not shown here). In this case, the other two types of crystal dislocations, with Burgers 
vectors 1/2 [110] and 1/2 [1 10], would form the dislocation walls at the twin boundaries. 

VII. CONCLUSIONS 

The Landau-Ginzburg theory developed in this paper represents the first step in the 
formulation of a mesoscopic model for martensitic phase transformations mediated by de- 
fects. In this framework the presence of dislocations induces incompatibility between the 
components of the elastic strain tensor and, therefore, the Saint- Venant condition no longer 
applies in media with defects. Instead, the elastic strains are related by an incompatibility 
constraint and the "strength" of this incompatibility is proportional to the gradients of the 
components of the Nye dislocation density tensor 12 . This incompatibility can be written in 
terms of densities of crystal dislocations in individual discrete slip systems which provides 
a clear "recipe" for coarse-graining microscopic information to the mesoscopic description. 
The incompatibility field is then completely determined and we have shown that its exis- 
tence leads to such interesting phenomena as the dependence of the order parameter texture 
on the density of dislocations or formation of correlated dislocations walls along the twin 
boundaries below T c . 

The presence of dislocations, i.e. the existence of a finite incompatibility field, introduces 
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new contributions to the free energy functional, in particular coupling of the order parameter 
with the incompatibility field and with the applied stress tensor. We have shown that 
the free energy minimum cannot be reached by minimizing independently all terms that 
depend on the order parameter. During this minimization the coupling between the order 
parameter and the incompatibility field introduces competition between these terms and the 
minimization of the free energy is thus inherently frustrated by the finite dislocation density. 
The order parameter field that minimizes the free energy subject to a given distribution of 
dislocations (i.e. strain incompatibilities) can be directly used to calculate internal strain 
and stress fields. For an arbitrarily anisotropic material, one can thus obtain the Peach- 
Koehler forces on the dislocations just by minimizing the free energy. These forces have 
been employed in the evolution equations for the density of individual variants of crystal 
dislocations (n s+ and n s ~) which takes the form of the Fokker-Planck equation. Hence, the 
conservation of the total Burgers vector is automatically satisfied. This procedure represents 
a simple self-consistent scheme in which the order parameter field is first calculated by 
minimizing the free energy subject to a given distribution of strain incompatibilities, i.e. 
densities of crystal dislocations. The corresponding internal stress field and the Peach- 
Koehler forces are then used to update the dislocation density using the discretized version of 
the Fokker-Planck equation. The new distribution of strain incompatibilities that correspond 
to the updated dislocation density is used again to minimize the free energy and this in turn 
provides the new order parameter field. 

To demonstrate the main features of the model developed in this paper, we considered 
two case studies that represent the limits of the theory with regard to the mobility of 
the dislocations. In the first case the four crystal dislocations with the 1/2(110) Burgers 
vectors were considered as immobile and, therefore, the free energy was subject to a fixed 
distribution of strain incompatibilities and minimized purely by the order parameter field 
62- We have shown that the nucleation of martensite and its spatial distribution depends on 
the density of dislocations. At low dislocation densities the minimization of the free energy 
is only weakly affected and the corresponding microstructure below T c is represented by 
martensite twins along the (110) directions. In contrast, at high dislocation densities the 
distribution of strain incompatibilities strongly affects the minimization of the free energy 
and eliminates the twinned microstructure. In the second case study we considered that all 
dislocations are completely mobile. We have demonstrated that below T c the free energy is 
minimized by forming alternating net dislocation walls at the mesoscale that decorate the 
twin boundaries between different variants of the martensite. The results are qualitatively 
consistent with the observations of correlated dislocation domains along twin boundaries in 
Ag 39 , Ni-Ti 7 and Fe-Pd 40 . 

Our objective here was mainly to demonstrate how continuum theory of dislocations 26 
can be incorporated into the mesoscopic free energy functional for displacive phase 
transformations' 1 which inevitably led to a number of simplifications. Many crystallographic 
details are beyond the resolution of the model and thus are included only in a coarse-grained 
manner. In particular, short-range interactions between crystal dislocations are replaced by 
continuous dislocation density fields n s+ and n s ~ that correspond to the two directions of 
the Burgers vector in the slip system s. The individual slip systems are assumed to be inde- 
pendent of each other and no reactions between dislocations moving in different slip planes 
are incorporated. No upper limit on the density of dislocations due to a finite minimum 
distance between dislocations of the same type is imposed at this point. Since no external 
loading is applied here, we do not consider the existence and operation of dislocation sources. 



19 



In the future, the model developed in this paper will be advanced to study the effect of 
dislocations on strain hardening and hysteresis, both of which play crucial roles in shape 
memory alloys. It serves as the basis for analytical solutions of elastic fringing fields at 
habit planes in the presence of dislocations and thus allows for a study of the role of defects 
on the size dependence of the twinning width in martensites. The concept of the defect- 
induced incompatibility of elastic strains is general and the approach pioneered in this paper 
may be applied to other phase transitions that are mediated by defects. Examples include 
the study of strain-induced polarization in ferroelectrics, strain-induced magnetization in 
ferromagnets or even the effect of long-range strain fields of dislocations on the recently 
discovered ferrotoroidal ordering 1(1 in ferrotoroidic materials. 
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APPENDIX A: TENSORIAL REPRESENTATIONS 



For convenience, following are tensorial representations of the operations that are used 
throughout this paper 1 '. Here, is the antisymmetric Levi-Civita tensor and A a tensor 
of rank two. 



divergence: 
curl: 

symmetric curl: 
incompatibility: 



V- A = A 



sym (V x A) = lie 



V x A = €i 
sym (V x j 

V x V x A 



^irs^jmn-^-sn,mr 



APPENDIX B: K-SPACE KERNELS 

In the following we write explicitly all the k-space kernels that are used to calculate 
the free energy and to perform the relaxational dynamics for the order parameter ei- The 
denominators in these kernels are identical and we evaluate them separately as 

d{k) = -j-fc 4 + -^k 2 x k 2 v . (Bl) 
v ' A x A 3 x y v ' 

The following kernels are used to calculate the secondary order parameter fields e± and e 3 
from the known primary order parameter field 62, the incompatibility field 7733, and the 
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external stress field ov,-: 



Qi(k) 



Si(fc) 



1 fc 4 A; 4 
1^ d{k) 

Ri(k) = 
1 fc 4 



, Q 3 (fe) 
l fc 2 



8 k x ky(k^, ky) 



A 1 d(k) 



A 3 d(k) 
Rz(k) = - 



8 k x ky 



Ti(k) 



A 1 d(k) 



8 Jv X fcjyfcj 



S 3 (k) 



n(k) 



i 



A, d(k) 
A 1 A 3 d(k) 



;b2) 



a 3 d(fe) 



- 1 



^3 d(fe) 

The nonlocal part of the free energy density, fnonioc, and its partial derivative, 
dfnonioc/de 2 , are expressed using 



An(k) 



" 2 



Bn(k) 



h 2 - h 2 



Cn(k) 



1 



d(k) ' lov ' d(k) ' lov ' d(k) 
and the kernels that depend on the above and the applied stress tensor are 

E AQ (k) = AiQi(fc)mi(fc) - A 3 Q 3 (k)m 3 (k) , 
E AR (k) = A 1 R 1 (k)m l (k) - A 3 R 3 (k)m 3 (k) , 
E A (fc) = ^[m^fc^ + Aalm^fc)] 2 , 

where: 

mi(fe) = S , i(fe)o- ii (fc)-T 1 (fc)t7 12 (fc) , 
m 3 (fc) = S^tydjjik) - T 3 (k)a 12 (k) , 



(B3) 



(B4) 



(B5) 



Here, the subscripts of each E imply which combination of the coefficients A+, Qi, Ri is used, 
and (Tjj = an + a 22 . 

Finally, the strain energy density, fi oa d, is expressed using the kernels that depend only 
on the applied stress tensor, a^, as: 

W Q {k) = \. {[Qi(fc) + 1] * n (k) + [Qxik) - 1] a 22 (k)} + Q 3 (k)a 12 (k) , 



(B6) 



W R (k) = -^R^kfajik) + R 3 {k)a 12 {k) , 
W s (fc) = —j=<Tjj(k)mi(k) - (Ji 2 (k)m 3 (k) . 



APPENDIX C: TEMPERATURE DEPENDENCE OF THE ELASTIC CON- 
STANTS FOR FE-30AT.%PD 

The temperature dependence of the elastic constants of Fe-30at.%Pd alloy, measured by 
Muto et al. and shown in Fig. 9, is approximated by linear relations given in Table I. The 
scaling laws are obtained by fitting the experimental data for Cn, (Cn + C\ 2 + 2C*44)/2 
and (Cn — C\ 2 )/2 that are linear within the range of temperatures considered in these 
measurements. The remaining parameters B, C, and K 2 are taken from Ref. 14 and are 
assumed to be independent of temperature: B = — 1.7 x 10 4 GPa, C = 3 x 10 7 GPa, 
K 2 /a 2 = 25 GPa. 
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TABLE I: Temperature scaling of the elastic constants and the coefficients used in the free energy 
(14). Units: T [K], Cy [GPa], A { [GPa]. 

C U (T) C 12 (T) C AA (T) 

104.82 + 0.1279T 229.09 - 0.3321T 71.73 + 0.0010T 



Ai(T) MT) A 3 (T) 

333.91 - 0.2042T -124.27 + 0.4600T 282.00 + 0.0344T 
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FIG. 9: Temperature scaling of elastic constants for Fe-30at.%Pd from Table I compared with the 
experimental data of Muto et al. . 
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49 The term "net dislocation" is used here in a loose sense and should not be thought of as a 
crystal dislocation. The reason is that the net Burgers vector is a vector sum of the Burgers 
vectors of individual crystallographic dislocations and, therefore, it can assume in general any 
orientation and magnitude. 

50 In the following, we will show that both elastic and plastic parts of the distortion tensor contain 
signatures of individual dislocations. To simplify the notation, we avoid the superscript e that 
often labels the elastic part of the distortion tensor. 

51 Tensorial representations of divergence, curl and incompatibility are given in Appendix A. 

52 The Saint- Venant law is satisfied for all irrotational strain fields, for example those due to point 
defects. 

53 There seem to be two equivalent definitions of this tensor in the literature. In the original 
Kroner's formulation, a = — V x (3 and the incompatibility of strains was then defined as 
r\ = V x e x V. Here, we adopt the convention of El-Azab whereby a = V x /3, and this 
allows us to write the incompatibility constraint in a more intuitive way as// = VxVxe. 

54 A similar approach is valid also in higher dimensions, where typically more than one primary 
order parameter is needed to identify the phase transition. 

55 In the isotropic case, C\\— C\ 2 = 2C44 and, therefore, only two elastic constants are independent. 
Hence, only two of the three coefficients A\, A 2 , ^3 in the free energy are independent since 
A 2 = A3/2. 

56 If we consider periodic boundary conditions, the integration by parts transfers the derivatives 
of the variations of e± and e^, to the derivatives of A, e.g. AV 2 <5ei becomes (V 2 A)<5ei, etc. 

57 If the incompatibility field 7733 vanishes, i.e. the medium is dislocation-free, and no external 
stress is applied, only the first term in (18) remains. In this case, one obtains 11 that e\ and e% 
can be expressed as functionals of the primary order parameter e 2 only. 

58 It is always possible to augment this deterministic relaxational dynamics by a stochastic term 
that is often characterized as white noise with zero mean and variance 2TksT. However, it 
can be shown that this noise plays an important role only at the temperatures slightly below 
T c , where thermal fluctuations may overcome the energy barrier between the two variants of 
the martensite and, therefore, cause switching between these variants. In this case, the twin 
boundary between individual martensite variants would not be sharp but rather diffuse. At low 
temperatures, i.e. well below T c , this barrier is large and the weak thermal fluctuations cannot 
cause this switching. Similarly, above T c , the free energy has one minimum corresponding to 
e 2 = and the thermal fluctuations would merely cause broadening of the distribution of e 2 . 
The goal of this paper is to give a proof of the principle and, for simplicity, the effect of thermal 
noise is not included. 
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